Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

High colored dissolved organic matter (CDOM) absorption in surface waters of the central-eastern Arctic Ocean: Implications for biogeochemistry and ocean color algorithms

  • Rafael Gonçalves-Araujo ,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

    rafael.goncalves.araujo@awi.de, rafaelgoncalvesaraujo@gmail.com

    Current address: Oceanographic Institute, University of São Paulo, Praça do Oceanográfico 191 –Cidade Universitária, Butantã, São Paulo, Brazil

    Affiliations Phytooptics Group, Physical Oceanography of Polar Seas, Climate Sciences Division, Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research, Bremerhaven, Germany, Faculty of Biology and Chemistry (FB-2), University of Bremen, Bremen, Germany

  • Benjamin Rabe,

    Roles Funding acquisition, Resources, Writing – review & editing

    Affiliation Physical Oceanography of Polar Seas, Climate Sciences Division, Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research, Bremerhaven, Germany

  • Ilka Peeken,

    Roles Data curation, Formal analysis, Funding acquisition, Methodology, Writing – review & editing

    Affiliation Polar Biological Oceanography, Biosciences Division, Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research, Am Handelshafen 12, Bremerhaven, Germany

  • Astrid Bracher

    Roles Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Supervision, Visualization, Writing – review & editing

    Affiliations Phytooptics Group, Physical Oceanography of Polar Seas, Climate Sciences Division, Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research, Bremerhaven, Germany, Institute of Environmental Physics, University of Bremen, Bremen, Germany

Abstract

As consequences of global warming sea-ice shrinking, permafrost thawing and changes in fresh water and terrestrial material export have already been reported in the Arctic environment. These processes impact light penetration and primary production. To reach a better understanding of the current status and to provide accurate forecasts Arctic biogeochemical and physical parameters need to be extensively monitored. In this sense, bio-optical properties are useful to be measured due to the applicability of optical instrumentation to autonomous platforms, including satellites. This study characterizes the non-water absorbers and their coupling to hydrographic conditions in the poorly sampled surface waters of the central and eastern Arctic Ocean. Over the entire sampled area colored dissolved organic matter (CDOM) dominates the light absorption in surface waters. The distribution of CDOM, phytoplankton and non-algal particles absorption reproduces the hydrographic variability in this region of the Arctic Ocean which suggests a subdivision into five major bio-optical provinces: Laptev Sea Shelf, Laptev Sea, Central Arctic/Transpolar Drift, Beaufort Gyre and Eurasian/Nansen Basin. Evaluating ocean color algorithms commonly applied in the Arctic Ocean shows that global and regionally tuned empirical algorithms provide poor chlorophyll-a (Chl-a) estimates. The semi-analytical algorithms Generalized Inherent Optical Property model (GIOP) and Garver-Siegel-Maritorena (GSM), on the other hand, provide robust estimates of Chl-a and absorption of colored matter. Applying GSM with modifications proposed for the western Arctic Ocean produced reliable information on the absorption by colored matter, and specifically by CDOM. These findings highlight that only semi-analytical ocean color algorithms are able to identify with low uncertainty the distribution of the different optical water constituents in these high CDOM absorbing waters. In addition, a clustering of the Arctic Ocean into bio-optical provinces will help to develop and then select province-specific ocean color algorithms.

1. Introduction

The Arctic Ocean basin receives 11% of the global freshwater input with its volume representing only 1% of the global ocean [1]. It obtains the largest amount of freshwater relative to its volume and therefore is the ocean most influenced by the continents on Earth. Together with the fresh water, high loads of terrestrial material (organic and inorganic; dissolved, colloidal and particulate) are introduced in that basin, in particular through the wide Siberian continental shelves [26]. By this the Arctic Ocean presents a large carbon reservoir and plays an important role in the planet’s carbon cycle. Besides, the Arctic environment has been experiencing the effects of ongoing global warming regarding permafrost thaw [7], changes in fresh water export [8,9] and decline of sea-ice extent [10,11] and volume [12]. The permanent loss of sea-ice may lead to an increase in light penetration in the Arctic surface layer [13] and to changes in the composition of phytoplankton assemblages [14], the overall primary production in the Arctic Ocean [15,16], and the degradation of terrestrial material transported to that basin [17,18].

Recent studies have pointed out regional differences in the Arctic Ocean with respect to biogeochemical parameters. For instance, shelf and open Arctic seas have shown to diverge in regards to the fluxes of biogenic matter [4] and export of terrigenous material [3]. Furthermore, while varying between shelf and open water in the Arctic [3,19], colored and fluorescent dissolved organic matter (CDOM and FDOM, respectively) content also differ between the western and eastern Arctic seas [6,2023]. Similarly, geographical differences in primary [24,25] and net community production [26], as well as in phytoplankton (e.g., dinoflagellates) and protist distribution [27,28] in the various basins of the Central Arctic Ocean have been observed. Such biogeographic patterns are likely related to hydrographic and sea-ice conditions within the region [28,29], denoting a strong coupling of physical and biogeochemical processes within the surface layers of the Arctic Ocean.

With the aforementioned effects of global warming and its impacts on the Arctic environment, improved monitoring and understanding of the current situation and changes in biogeochemical parameters are necessary. The optical properties of dissolved organic matter are reliable water mass tracers in the Arctic Ocean according to reports based on in situ [22,30] and remote sensing data [31] and has also been proven to be useful on monitoring small scale changes in coastal environments [32]. Biogeochemical parameters such as chlorophyll-a (Chl-a) and CDOM can be determined (and be estimated for primary production) in surface waters by ocean color remote sensing. Furthermore, the Arctic Ocean is a unique ocean where, even in pelagic waters, the non-water light absorption in the surface layer being dominated by CDOM [21] which does not co-vary with Chl-a. Opposed to that, the latter is assumed by empirical ocean color algorithms. Hence, these algorithms lead to an overestimation of Chl-a [3335] and overall poor performance in the Arctic [36]. Improvement of algorithms for the Arctic Ocean is challenging given the difficulties to sample for validation data in those waters, in particular, on the remote Siberian shelves [37]. Several studies have addressed the quality of the estimates from ocean color algorithms in the western Arctic ocean [19,21,3335,38,39]. Regionally tuned algorithms provided improved estimates related to global algorithms in the western Arctic [34,38]. In that same region, semi-analytical algorithms obtained even better estimates of Chl-a [21,33]. Besides Chl-a, semi-analytical algorithms can also retrieve CDOM in that region with low uncertainty [19,21,33].

Here, we focus on optical and hydrographic sampling in the central-Eastern Arctic, an area up to now hardly evaluated [39] for the application of satellite and in situ optical measurements to monitor the surface biogeochemistry of the Arctic Ocean. The objectives of this study are twofold: first, we aim to obtain a characterization of the non-water absorption constituents in the surface waters in the Central-Eastern Arctic. Those properties were tested whether they reproduce hydrographic and geographic patterns (or units). As a second objective, we evaluate empirical and semi-analytical ocean color algorithms commonly applied to studies in the Arctic Ocean and compare their performances. Given the novelty of the results presented in this study, it contributes to the growing Arctic remote sensing research, which has been so far mostly devoted to the western Arctic Ocean. Moreover, as already pointed out [37], the sampling effort for the Arctic Ocean is still very low compared to other ocean basins and more studies are required to improve the ocean color estimates for that basin. Finally, it is important to stress that whilst ocean color sensors are not able to monitor under very low (or no) illumination and cloudy conditions and ice-covered regions, in situ bio-optical measurements in those regions are crucial for improving biogeochemical models; however, such measurements are very scarce in the central and eastern Arctic Ocean. Furthermore, results on in situ bio-optical and biogeochemical properties are important for calibrating sensors coupled to autonomous platforms (e.g. satellites; gliders; Autonomous Underwater Vehicles, AUVs; Ice-Tethered Platforms, ITPs; etc.). In the future, those sensors will measure in situ biogeochemical properties enabling monitoring on high spatial and temporal resolution and coverage in the Arctic Ocean [e.g., ITPs [4042]].

2. Methods

2.1. Sampling

The ARK XXVI-3 (PS-78) cruise was conducted in shelf and open waters through the central-eastern Arctic Ocean from 5th August to 6th October 2011 onboard the R/V Polarstern. Temperature and salinity profiles were acquired with a CTD attached to a rosette system at 110 oceanographic stations [43] (Fig 1A). Surface water samples for analysis of dissolved organic matter, particulate matter and chlorophyll-a (Chl-a) were taken using Niskin bottles attached to the rosette system at 62 stations (Fig 1B). As observed in Fig 1C and 1D, most of the sampled area was covered by sea-ice. No specific permissions were required for these locations/activities given that sampling was performed out of the 200 Mile zone. Data are available at https://doi.org/10.1594/PANGAEA.867532.

thumbnail
Fig 1. Study region and sea-ice conditions.

ODV maps [44] showing the sampling stations occupied during the ARK-XXVI/3 (PS-78) cruise where CTD casts (a), water sampling and hyperspectral radiometric measurements (b) were performed. Arrows in (a) represent the main surface circulation patterns in the Arctic Ocean colored as follows: major rivers (green); inflowing currents (red); out flowing currents (blue) [45]. AMRSR-2 sea-ice concentration (http://meereisportal.de) for August (c) and September (d) 2011.

https://doi.org/10.1371/journal.pone.0190838.g001

2.2. Particulate absorption analysis

Water samples for particulate absorption analysis were filtered on GF/F filters (0.7 μm pore size), shock-frozen in liquid nitrogen and stored at –80°C until laboratory analysis at the Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research. Measurements were carried out on a dual-beam UV/VIS spectrophotometer (Cary 4000, Varian Inc.) equipped with a 150 mm integrating sphere (external DRA-900, Varian, Inc. and Labsphere Inc., made from SpectralonTM) using a quantitative filterpad technique [46]. The filters were placed inside and at the center of the integrating sphere using a center-mount filter holder perpendicular to the light beam. A wavelength scan from 300 to 850 nm with a resolution of 1 nm (slit width 2 nm, scan rate 150 nm min1) was performed, when the reflectance ports were covered with SpectralonTM reflectance standards. The baseline was recorded beforehand with a clean, dry filter, and a filter, which was soaked for more than 30 min in freshly produced Milli-Q water, was taken as a reference. The absorption coefficient was calculated from the transmittance [T(λ)], which is derived from the optical density (OD) measurements, using a path length amplification factor of 4.5 (β = 1/4.5) [47] following the equation: (1) (2) where V is the filtrated sample volume in m3 and A the filter clearance area in m2. Results from the original filter gave total particulate absorption, ap. The algal pigments were bleached with NaOCl [48,49] to determine the absorption by detrital material (or non-algal particles), hereafter referred to as non-algal particles (aNAP). The bleached filters were measured following the procedure described above. The particulate absorption of phytoplankton at each wavelength (λ) [aph(λ)] was obtained by subtracting aNAP from ap.

2.3. Dissolved organic matter absorption analysis

Water samples for DOM analysis were filtered through prerinsed 0.2 μm filters immediately after sampling and stored in amber glass vials in dark at 4°C until analysis in laboratory at the Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research after the cruise. CDOM was analyzed with an Aqualog® fluorescence spectrometer (HORIBA Jobin Yvon, Germany) using freshly produced Milli-Q water as reference. CDOM absorbance spectra measurements (260–600 nm) were blank-corrected and a baseline correction was applied at 600 nm, assuming negligible CDOM absorption at that wavelength. CDOM absorbance was further converted into Napierian absorption coefficient [aCDOM(λ)], obtained from the given equation: (3) where A(λ) is the absorbance at specific wavelength and L is the cuvette path length in meters. a is generally adopted as a proxy for assessing the CDOM content in a given water sample and in this study it is presented in the visible [440 nm—aCDOM(440)] and UV [350 nm—aCDOM(350)] bands. aCDOM(440) was chosen given its application to ocean color remote sensing [50,51] and to make it comparable with the particulate matter absorption coefficients [52]. The UV band aCDOM(350) was determined in this study due to its correlations to DOC and lignin concentrations and to permit comparison with earlier results [6,20,53].

2.4. Chlorophyll-a analysis

For measuring the photosynthetic pigment Chl-a one liter of seawater samples were taken from Niskin bottles, immediately filtered on GF/F filters, frozen in liquid nitrogen, and stored at -80°C until further analyses by high-performance liquid chromatography (HPLC) at the home laboratory of the Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research after the cruise. The samples were measured using a Waters 600 controller equipped with an auto sampler (717 plus), a photodiode array detector (2996) and the EMPOWER software. Chl-a was analyzed by reverse-phase HPLC using a VARIAN Microsorb-MV3 C8 column (4.6 3 100 mm) and HPLC-grade solvents (Merck). The solvents gradient and routine of analysis are fully described in Taylor et al. [54]. Chl-a concentrations were quantified based on peak area of the external standard, which was spectrophotometrically calibrated using extinction coefficients published in Jeffrey et al. [55].

2.5. Radiometric measurements

Underwater optical light fields were assessed through radiance and irradiance profiles obtained with hyperspectral radiometers (RAMSES, ARC-VIS and ACC-VIS, respectively, TriOS GmbH, Germany). The instruments cover a wavelength range of 320 nm to 950 nm with an optical resolution of 3.3 nm and a spectral accuracy of 0.3 nm. Measurements were collected with sensor-specific automatically adjusted integration times (between 4 ms and 8 s). 16 radiometric profiles (Fig 1B) were collected simultaneously with the CTD profiles down to a maximum depth of 100 m. At each profile, measurements of upwelling radiance (Lu) and downwelling irradiance (Ed) were performed. One of the in-water sensors was equipped with inclination and pressure sensors. To avoid ship shadow, the ship was oriented such that the sun was illuminating the side where the measurements have taken place.

The radiometric measurements were performed out of the ship’s shadow and during clear sky or nearly clear sky conditions; this was checked based on the ship’s global radiation sensor data ensuring low variation of the incoming sunlight. For the in-water data, the inclination in either dimensions was smaller than 14° [35]. During the acquisition of the profiles, stops (varying from 30 to 60 s) were performed within a 10 m depth interval. These data were then averaged in discrete intervals of 5 and 10 m for 0–30 m and below, respectively, and were further processed following the NASA protocols [56]. As surface waves strongly affect measurements in the upper few meters, deeper measurements that are more reliable to be used can be further extrapolated to the sea surface [56]. Analogously to Stramski et al. [57] a depth interval was defined (z’ = 10 to 30 m) to calculate the vertical attenuation coefficients for downwelling irradiance and upwelling radiance, [i.e. Kd(λ,z’) and Ku(λ,z’), respectively]. With Kd(λ,z’) and Ku(λ,z’), the subsurface irradiance Ed(λ, 0 m) and radiance Lu(λ,0 m) were extrapolated from the profiles of Ed(λ,z) and Lu(λ,z).

For the calculation of the remote sensing reflectance [Rrs(λ)], the subsurface Lu(λ, 0 m) and Ed(λ, 0 m) were propagated through the water-air interface by applying a transfer coefficient of 0.519 [57]. Rrs(λ) was then calculated: (4)

2.6. Ocean color algorithms

In this study we evaluated the performance of ocean color algorithms to derive Chl-a, aCDOM(λ), adg(λ) [the sum of aCDOM(λ) and aNAP(λ)] and aph(λ). Firstly, we tested different empirical algorithms, which are used to derive Chl-a from band ratios of Rrs. These algorithms are frequently applied to the Arctic Ocean. Here their Chl-a retrievals were obtained using Rrs from the 16 stations as input and then compared with in situ measured Chl-a. The MODIS OC3M is a global algorithm, which is determined as a function of three Rrs band ratios [58]. The global SeaWiFS OC4V6 [58,59] and the regional Arctic OC4L [38] algorithms, nevertheless, use a four-band ratio approach. These algorithms are expressed as follows: (5) (6) (7) where R is the base 10 logarithm of the maximum band ratio, whichever is the greatest of Rrs(443)/Rrs(555), Rrs(490)/Rrs(555), and Rrs(510)/Rrs(555); R’ is the same as R but it considers the greater of the two band ratios Rrs(443)/Rrs(551) and Rrs(488)/Rrs(551); and the coefficients a, b, c, d, e, a1, b1, c1, d1, e1, a2, and b2 are empirically derived values. Additionally, the performance of modifications to the global OC3M and OC4V6 algorithms developed for the western Arctic Ocean [34] hereafter OC3M-mod and OC4V6-mod, respectively, was evaluated. The coefficients for those regional algorithms are given below:

  • OC3M-mod: a3 = –0.32, b3 = –2.33, c3 = 4.02, d3 = –31.64, e3 = 48.54;
  • OC4V6-mod: a4 = –0.35, b4 = –1.52, c4 = –2.44, d4 = –12.80, e4 = 30.48.

Apart from the empirical ocean color algorithms, two semi-analytical algorithms (SAA) were tested. First, we used the Generalized Inherent Optical Property model (GIOP) [60,61], for simplicity further named GIOP, using settings applied for the western Arctic [33] to allow comparison with the results from that study. In short, GIOP is a spectral matching inversion model, which applies non-linear least square methods to retrieve three eigenvalues [aph(443), adg(443) and the particles spectral backscattering coefficient–bbp(555)]. GIOP can also estimate Chl-a from aph(443), by using the factor 0.055. As in Chaves et al. [33], we used the GIOP applied to in situ Rrs(λ) at the SeaWiFS/MODIS-Aqua operational wavelengths (412, 443, 490, 510, 555 and 670 nm). Besides the GIOP, a modification of the Garver-Siegel-Maritorena (GSM) SAA [62,63] for retrieving aCDOM(λ) in the Arctic Ocean [19] was used. This algorithm was developed based on a parametrization of absorption properties using data from the western Arctic. In short, it enables the separation of aNAP(λ), and therefore aCDOM(λ), from adg(λ) by applying the parametrization of aNAP(λ) related to the particle backscatter at 555 nm [bbp(555)] [35].

To summarize, in this study we evaluate the following retrievals from ocean color algorithms:

  • Chl-aOC3M [58];
  • Chl-aOC4V6 [58,59];
  • Chl-aOC3M-mod and Chl-aOC4V6-mod [34];
  • Chl-aOC4L [38];
  • aph(λ)GIOP, adg(λ)GIOP, and Chl-aGIOP [60,61];
  • adg(λ)Mat and aCDOM(λ)Mat [19].

2.7. Statistical analysis

Hierarchical cluster analysis using simple average linkage and Euclidean distance method was applied to classify both, a matrix containing hydrographic and inherent optical properties (IOPs) bulk properties (hereafter environmental matrix) and a matrix consisting of hyperspectral apparent optical properties (AOPs), into hydrographically (and geographically) coherent groups [64]. The purpose of applying such an approach was to test whether the AOPs, or a combination of hydrography and IOPs, are capable to trace the origin of waters masses as previously suggested in other regions of the Arctic Ocean [22,31]. The first matrix, named environmental matrix, consisted of surface measurements of temperature, salinity, aCDOM(443), aNAP(443) and aph(443), which were normalized prior to analysis, by subtracting the mean value and then dividing by the standard deviation. The environmental matrix consisted of parameters, which in future can be derived from sensors mounted together on autonomous platforms. The hyperspectral AOP matrix consisted of the second derivative of Rrs(λ)/Rrs(555) that computes the changes in curvature of a given spectrum over a sampling interval or band separation [65]. For that, Rrs spectra were interpolated to the optimal range for band separation (435–510 nm), and a smoothing filter window of 27 nm was used for the derivative calculations [65]. Given that bio-optical parameters are generally log-normally distributed in natural environments [66] and also in this study, power functions were applied to evaluate the correlation between pairs of bio-optical parameters [67,68].

Kruskal-Wallis H tests were applied to compare variables between pairs of clusters, after being normality-tested with the Kolmogorov-Smirnov test. To evaluate the performance of the ocean color algorithms, r2, slope, intercept, root mean square error (RMSE), mean absolute error (MAE) and bias for each pair of variables were determined. The RMSE, MAE and bias were calculated as follows: (8) (9) (10) where Y is the retrieved parameter (e.g. Chl-a, adg, aph) and X is the correspondent in situ measured parameter.

3. Results and discussion

3.1. Hydrography

Based on the temperature and salinity profiles five water masses were identified within the surface layer (0−200 m) of the sampled area, which are in agreement with previous studies in the region [45,69]: Upper Halocline Water (UHW), Barents Sea Branch Water (BSBW) and Laptev Sea Shelf Water (LSSW) at the surface; and Lower Halocline Water (LHW) and Atlantic Water (AW) in the beneath layer (Fig 2A). Surface waters of the central Arctic were occupied by the UHW, whereas BSBW and LSSW were observed at surface in pelagic and shelf waters in the Laptev Sea, respectively. Most of the Arctic surface waters are of Atlantic origin and are progressively modified at higher latitudes by heat exchange with the atmosphere, river runoff, melt water in summer and salt rejection in winter [70]. The UHW was observed within the upper polar mixed layer (PML, ~40 m) and was characterized by subzero temperatures and a thin upper mixed layer (10–20 m) due to dilution from sea-ice melt. This water mass can be further divided into two origins, given differences in the salinity ranges [45]: the UHW from the Canadian Basin with the Beaufort Gyre (cUHW, with salinity <32.5), and the UHW from the Eurasian Basin (eUHW, with salinity between 32 and 34). The BSBW observed in the surface waters of the Laptev Sea is transported from the Kara Sea through the Vilkitsky Strait by the Vilkitsky Strait Current [69,71,72]. High temperatures and low salinity characterize the LSSW which is strongly influenced by the Lena River outflow [73]. The highest temperature (3.67°C) and lowest salinity (23.8) values were observed within LSSW. The lowest temperature (−1.86°C), on the other hand, was associated with LHW whereas the highest salinity (34.9) was observed within AW located in the deepest sampled layer, generally below 100 m (Fig 2B and 2C).

thumbnail
Fig 2. Hydrography in the surface central and eastern Arctic Ocean.

(a) T-S diagram with depth (m) as color bar. Surface distribution of temperature (°C) (b) and salinity (c) with the approximate occupation of the water masses with the PML within the study region. Produced with Ocean Data View [44].

https://doi.org/10.1371/journal.pone.0190838.g002

3.2. Absorption coefficients of water constituents

Phytoplankton absorption coefficients [aph(443)] were highly correlated with the absorption coefficients of NAP [aNAP(443)] (r2 = 0.95; p<0.0001; n = 62). aph(443) ranged from 0.01 to 0.06 m-1 whereas the aNAP(443) varied between 0.0004 and 0.04 m-1. The highest aNAP(443) values were associated with sites close to the shelf break, denoting the continent as its main source, reaching its maximum within the LSSW, in similar ranges as previously reported [21]. Relatively high values of aph(443) were observed close to the Laptev Sea shelf break, as for aNAP(443), however, the highest aph(443) values were obtained for the Nansen Basin. Further discussion on the spatial variability of those parameters is presented in Section 3.3.

To investigate the correlation of aph(443) and ap(443) with Chl-a a power function was applied [67,68]. Both aph(443) (Fig 3D) and ap(443) were highly correlated to Chl-a, however, as expected, the correlations for aph(443) were higher. The power functions for aph(443) and ap(443) in relation to Chl-a concentration obtained in this study are given below: (11) (12)

thumbnail
Fig 3. Absorption coefficients in the surface waters of central and eastern Arctic Ocean.

(a) Surface distribution of aCDOM(443) (m-1) and correlation between aCDOM(443) and aCDOM(350) (inset graph); produced with Ocean Data View [44]. (b) Ternary plot denoting the contribution of the non-water absorbers [aCDOM(443), aph(443), aNAP(443)] to total non-water absorption [atw(443)] at surface; color bar indicates salinity. (c) Station 207 (indicated by the arrow in a) as example of atw(λ), aCDOM(λ), aph(λ) and aNAP(λ) spectra (m-1). Dashed line indicates the position of 443 nm. (d) Correlation between Chl-a (mg m-3) and aph(443) (m-1); for the colors, please refer to Fig 4.

https://doi.org/10.1371/journal.pone.0190838.g003

As shown in Fig 3D, the correlation between aph(443) and Chl-a was comparable to the one found by Bricaud et al. [68], and to other expeditions conducted in the Fram Strait and adjoining seas [7476]. The consistency between these results thus, reiterates the applicability of such absorption measurements in the VIS-range as a proxy to retrieve Chl-a concentrations. The correlation between aph(443) and Chl-a observed in this study presents, in turn, a slight deviation from the trends reported for the Canadian Basin [35]. Since the authors claim that their deviation from the global average is caused by a different pigment packaging effect and/or pigment composition, in their specific region: this would in turn mean that for our data set the phytoplankton composition and light adaptations follow more the global average.

CDOM absorption coefficients in the visible and UV wavelength ranges [aCDOM(443) and aCDOM(350), respectively] were highly correlated (r2 = 0.99, p<0.0001) and ranged from 0.02 and 0.19 m-1 to 1.14 and 4.42 m-1, respectively (Fig 3A). The highest aCDOM(443) values [aCDOM(443)>1 m-1] were observed in the Laptev Sea associated to the LSSW, with values in similar ranges as previously reported for those waters [21,50,53]. High aCDOM values [~0.5 m-1 for aCDOM(443)] were observed in the central Arctic, which have been shown to have a high terrestrial signal [77], likely associated to transport of high-DOM Siberian Shelf waters [2] within the Transpolar Drift. That high-DOM signal can be traced even after significant removal during the transport of those waters to the central Arctic [20] and in the Fram Strait [22,23,30,7476,78,79]. The lowest aCDOM values [<0.2 m-1 for aCDOM(443)] were observed in the Beaufort Gyre and Amundsen and Nansen basins which is related to the influence of waters from the Norwegian and Barents Sea [45] that have a very low DOM content [78,80]. The low aCDOM(443) observed in the Beaufort Gyre corroborates the well-known DOM decrease towards the center of oligotrophic oceanic basins and gyres, where aCDOM(443) values tend to be close to zero [81,82]. Furthermore, aCDOM(443) differed significantly (p<0.001) between pelagic samples from the Central Arctic and Beaufort Gyre. Likewise, a recent study reported higher aCDOM(443) values in shelf waters of Eurasian basin in comparison to the Canadian basin [21]. Such a difference in the DOM background between the two basins is likely a reflection of the higher loads of DOM from Siberian Rivers [6,20]. Moreover, the differences between DOM from Eurasian and Canadian basins can be also detected in the intensity of visible DOM fluorescence, which can further distinguish the origins of fresh water exiting the Arctic Ocean [22].

We computed the total non-water coefficient absorption spectra [aCDM(λ)] as follows: (13) atw(λ) was strongly correlated with aCDOM(λ) (p<0.0001) in the UV and VIS (violet-blue, mostly) wavelength ranges, suggesting CDOM as the major absorber component of the surface waters in that spectral regions through the entire sampling area (Fig 3C). Such dominance of CDOM is also clear when looking at the relative proportion of the non-water absorbers [aph(443), aNAP(443) and aCDOM(443)] to atw(443) (Fig 3B), which shows that all sampled stations are classified as CDOM-dominated [52]. The proportion of aCDOM(443) was high (Table 1), with it contributing to over 50% at all sampled stations, reaching a maximum contribution of 99% to atw(443). Similar averaged values (0.85 ± 0.07) for the aCDOM(443) contribution to atw(443) was observed in a recent study conducted in the Eurasian Basin [21]. Nearly as high contributions of CDOM were also reported for the Canadian Basin (Chukchi Sea: 0.74±0.14; Western Arctic: 0.76±0.11) in that same study [21]. Our study shows that CDOM is not only the major non-water absorber in the western Arctic and shelf seas of the eastern Arctic [21,35], it also strongly dominates the non-water absorption in the central Arctic. Dominance of CDOM to the total non-water absorption has been primarily reported to coastal environments, classified as “Case-2 waters” because CDOM (and also NAP) does not covary with Chl-a [83,84]. Oligotrophic pelagic systems (as the Central Arctic), on the other hand, are generally characterized as “Case-1 waters”, where Chl-a is thought to be the dominant absorber and covaries with CDOM [83,84]. That assumption, however, is not applicable to the pelagic Arctic Ocean, whose non-water absorption is clearly dominated by CDOM that, in turn, does not covary with Chl-a (r2 = 0.01). Such an absence of covariance between CDOM and Chl-a has been already reported for the Labrador Sea [85] and Western Arctic Ocean [35,86]. The correlation observed in this study (r2 = 0.01) was, nevertheless, the weakest observed for the entire Arctic Ocean, and could be related to the greater contribution and variability of CDOM to the total non-water absorption in our investigated waters. Finally, aNAP(443) contribution to atw(443) was the lowest found for the Arctic waters, being likely negligible compared to aCDOM(443) contributions (Table 1).

thumbnail
Table 1. Relative absorption of non-water absorbers.

Averaged contribution of the absorption coefficients for each of the non-water absorbers (at 443 nm) to atw(443) in this and other studies carried out in different regions.

https://doi.org/10.1371/journal.pone.0190838.t001

3.3. Geographic clustering

Hierarchical cluster analysis was applied to the environmental matrix [temperature, salinity, aCDOM(443), aNAP(443) and aph(443)] to classify the sampling sites according to coherent groups with respect to hydrography and non-water absorption. A total of seven major clusters were identified and those were used to divide the study area into five distinct geographic zones (Fig 4): Laptev Sea Shelf, Laptev Sea (pelagic), Central Arctic/Transpolar Drift, Beaufort Gyre and Nansen Basin. Those zones were easily discriminated based on the surface values of the environmental matrix. The average and standard deviation of the analyzed parameters for each cluster are presented in Table 2. In short, cluster 1 characterized the surface Laptev Sea shelf waters, influenced by the Lena River outflow, with high temperature, low salinity, moderate aph(443) and the highest values of CDOM and NAP, in agreement with previous reports for that region [2,5,21,53,88]. Cluster 6 was composed by stations located in the pelagic and western domain of the Laptev Sea, with influence of shelf waters from the Kara Sea [71,72]. Those waters presented high temperatures, relatively low salinity and moderate values of aCDOM(443), aph(443) and aNAP(443). Clusters 2 and 5 united the stations located in the Central Arctic, over the Transpolar Drift stream [45], where the Arctic shelf waters with relatively low salinity and high aCDOM(443) are transported along the Arctic Basin [20]; however cluster 5 seems to be a transitional zone, with less influence of Arctic shelf waters, exhibiting lower aCDOM(443) and higher aph(443) compared to cluster 2. Cluster 3 grouped the stations located in the Beaufort Gyre. Those lower salinity waters [89] presented near freezing temperature and very low non-water absorption was observed, with aCDOM(443) and aph(443) exhibiting the lowest values among the seven clusters. These results corroborate previous findings showing Canadian Basin water with low Chl-a and primary production [16,25], as well as lower DOM content [6,2022], in comparison to the Eurasian Basin. Finally, the clusters 4 and 7 grouped the stations located in the Nansen and Amundsen basins, with influence of waters advected from the North Atlantic Ocean and Norwegian Sea. Those waters were characterized by the lowest temperatures, the highest salinity, low aCDOM(443) and aNAP(443), as reported for the waters of the Atlantic inflow to the Arctic in the Fram Strait [79]. On the other hand, aph(443) (and Chl-a) values within that cluster were the highest, likely explained by the advection of nutrient rich Atlantic water [90] that stimulates phytoplankton growth. Clusters 4 and 7 differed from each other only regarding the aph(443) (and Chl-a) values, with the highest values being observed in cluster 7. High aph(443) (and Chl-a) observed in the Nansen and Amundsen basins can be attributed to the high transmittance of light in those waters primarily due to the development of melt ponds in the sea-ice [13], which increases primary production in those areas [24].

thumbnail
Fig 4. Clustering based on environmental matrix.

(top) Dendogram (cophenetic correlation coefficient: c = 0.91) for sampling stations based on surface normalized values of an environmental matrix containing hydrographic and IOP parameters: temperature, salinity, aCDOM(443), aNAP(443) and aph(443). (bottom) ODV map [44] showing the position of each station according to the classification based on the hierarchical clustering. Inset graph shows the correlation between aCDOM(443) and salinity colored with respect to the clusters; black line indicates the best fit (p<0.01).

https://doi.org/10.1371/journal.pone.0190838.g004

thumbnail
Table 2. Average of parameters for the geographic clusters based on the environmental matrix.

Averaged values ± standard deviation of hydrographic/IOP parameters and geographic region for each of the clusters presented in Fig 4. Geographic regions acronyms: BG (Beaufort Gyre); EB (Eurasian Basin–Amundsen and Nansen basins); LS (Laptev Sea); LSS (Laptev Sea Shelf–Lena river influenced); TPD (Transpolar Drift).

https://doi.org/10.1371/journal.pone.0190838.t002

To test whether hyperspectral remote sensing information is capable of detecting hydrographic and bio-optical variability we have also applied hierarchical cluster analysis to hyperspectral Rrs [in this case, the 2nd derivative of Rrs(λ)/Rrs(555); see section 2.7]. Despite the low number of sampled stations (n = 16), the analysis yielded satisfactory results (cophenetic correlation coefficient: c = 0.87) and two main clusters were isolated (Fig 5). The partition based on hyperspectral data shows some similarities with the one provided by clustering the environmental matrix (see Fig 4). Cluster I comprised the Rrs spectra (i.e. stations) with lower aCDOM(443), located mainly in the Nansen and Amundsen basins and North Laptev Sea, under influence of waters from the North Atlantic, Norwegian Sea and also from Kara Sea. This cluster corresponds to the clusters 6 and 7 (and two stations of the transition cluster 5), with relatively low aCDOM(443) and influence of waters advected from the abovementioned regions. Additionally, the only station from the Beaufort Gyre, which also presented low aCDOM(443), was included in this same cluster I. Here we speculate that given the low number of stations performed, the multivariate analysis may not be able to solve such variability and grouped all the low CDOM spectra into one unique cluster. However, with an increased number of sampling stations, such variability would be easier to be detected in Rrs spectra. Cluster II isolated Rrs spectra from stations with high aCDOM(443) and lower Rrs (Fig 5 and Table 3), located in the central Arctic and close to the Laptev Sea shelf (Fig 5). Its corresponding environmental clusters are mainly the clusters 2 and 5, which were under influence of the shelf waters transported within the Transpolar Drift [45]. No stations of clusters 3 and 4 were sampled for hyperspectral remote sensing information.

thumbnail
Fig 5. Hyperspectral AOP clustering.

(a) Dendogram (cophenetic correlation coefficient: c = 0.87) for sampling stations based on hierarchical cluster analysis applied to the 2nd derivative of Rrs(λ)/Rrs(555) (following Torrecilla et al., 2011). (b) ODV map [44] showing the position of each station according to the classification based on the hierarchical clustering. (c) 2nd derivative of normalized hyperspectral remote sensing reflectance, Rrs(λ)/Rrs(555), with respect to the wavelength range of 435–510 nm (following Torrecilla et al. [65]). (d) Normalized hyperspectral remote sensing reflectance, Rrs(λ)/Rrs(555) in the visible wavelength range. Colored circles in (a) refer to the environmental clusters presented in Fig 4. Colors in (c) and (d) are in accordance with the clusters presented in (a) and (b).

https://doi.org/10.1371/journal.pone.0190838.g005

thumbnail
Table 3. Hydrographic and IOP parameters for the geographic clusters based on hyperspectral AOP measurements.

Averaged values ± standard deviation of geophysical parameters for each of the clusters presented in Fig 5.

https://doi.org/10.1371/journal.pone.0190838.t003

3.4. Arctic bio-optical provinces

The results provided by hierarchical cluster analyses in this study (see Figs 4 and 5 and Tables 2 and 3) show that hydrographic data and non-water absorption, but also hyperspectral AOPs (e.g. Rrs spectra) are applicable tools for characterizing surface waters (geographic zones) with differing surface biogeochemical properties, even in waters where non-water absorption is strongly dominated by CDOM, such as the Arctic Ocean [21]. Similarly, a recent study applied hierarchical cluster analysis to the spectral particulate backscattering-to-absorption ratio in the western Arctic allowing the partitioning of optically-distinct clusters of particles assemblages, which, in turn, reflect difference in the characteristics of particle concentration, composition, and phytoplankton taxonomic composition and size [91]. Furthermore, given the coupling between hydrographic and bio-optical properties, one can further suggest those clusters as bio-optical units or provinces. Bio-optical provinces based on HCA applied to IOPs and AOPs have shown to be reliable describers of Longhurst provinces [92] in the Atlantic Ocean [54]. On the other hand, almost the entire Arctic Ocean is classified as a unique ecological province, the Boreal Polar Province (BPLR), within the Polar Biome [92]. That same author suggested that there might be spatial variability between shelf and pelagic ecosystems (as well as in the marginal ice zones) within the BPLR, however it is very difficult to sustain an adequate description of smaller units, given the difficulty to access the northern seas. Along with that, differences among the Arctic Seas have been already reported, for instance, with respect to export of biogenic matter [4], number of dinoflagellates species [27], protist diversity [28,29], and primary production [24,25,93]. Those studies, therefore, reinforce the existence of distinct biogeographic units in the Arctic Ocean and further implementation of a biogeographic characterization in the region is of great importance to improve the current understanding about the Arctic environment. The determination of such biogeographic zones would guide future strategies for Arctic monitoring and ecosystem modeling, leading to a more accurate understanding of the ecosystem functioning and biogeochemical stocks, as well as on the prediction of future scenarios with regards to climate change. Finally, to build on that, based on the results presented by our quasi-synoptic sampling through the central-eastern Arctic Ocean, we therefore propose an overall classification of the sampling sites into five major bio-optical provinces. Those sites were defined based on a combination of hydrographic characteristics and IOPs, but also considering the outcome of the hyperspectral AOP matrix clustering. They are classified as follows (Fig 6):

  1. Laptev Sea Shelf: strongly influenced by the Lena River outflow, is primarily characterized by low salinity, high temperature and aNAP(443), moderate aph(443), Chl-a and very high aCDOM(443);
  2. Laptev Sea: low influence of Lena River outflow, however with contributions of waters advected from the Kara Sea; presents relatively low salinity, relatively high temperature and moderate levels of aCDOM(443), aph(443), Chl-a and aNAP(443);
  3. Central Arctic/Transpolar Drift: characterized by shelf waters transported within the Transpolar Drift, it has very low temperatures and relatively low salinity, aph(443) and Chl-a; however with high aCDOM(443) and very low aNAP(443);
  4. Beaufort Gyre: the waters with lower non-water absorption; they present low temperature and salinity, together with very low values of aCDOM(443), aph(443), Chl-a and aNAP(443);
  5. Eurasian/Nansen Basin: region influenced by waters advected from the Atlantic Ocean and Norwegian Sea, those waters present the highest salinity and near freezing temperature, with very low aCDOM(443) and aNAP(443), and the highest aph(443) and Chl-a levels due to high transmittance through sea-ice.
thumbnail
Fig 6. Arctic bio-optical provinces.

Distribution of the five Arctic bio-optical provinces defined in this study based on HCA applied to surface hydrographical, IOP bulk and hyperspectral AOP data.

https://doi.org/10.1371/journal.pone.0190838.g006

3.5. Evaluation of ocean color algorithms

Fig 7 shows the evaluation of the current global empirical ocean color algorithms OC3M and OC4V6 (and their regional adaptations for the western Arctic, OC3M-mod and OC4V6-mod) frequently applied to the Arctic Ocean, as well as the Arctic OC4L algorithm, which is designed to be applicable to high northern latitudes. When considering all sampled stations, the five empirical algorithms failed in retrieving Chl-a from Rrs bands, and a general inverse correlation with in situ Chl-a was observed (Fig 7 and Table 4). Furthermore, despite the relatively low RMSE observed for OC4V6-mod and OC3M-mod, all the band-ratio algorithms applied in this study appeared to attribute CDOM absorption to phytoplankton absorption (Fig 7D and Table 4). Such CDOM-biased retrievals from empirical Chl-a ocean color algorithms have already been reported for the western Arctic [3335]. This is attributed to the fact that CDOM is the greatest absorber at 443 nm over the entire sampled region (see Fig 3 and Table 1). As pointed out by Chaves et al. [33], excess aCDOM(λ)–that is assumed to co-vary with Chl-a–produces lower maximum band ratios [Rrs(443>490>510/555)], thus resulting in overestimation of Chl-a (see Fig 7C).

thumbnail
Fig 7. Evaluation of empirical ocean color algorithms frequently applied to the Arctic Ocean.

(a) Chl-a estimated by empirical algorithms (mg m-3; indicated by different colors) versus in situ Chl-a (mg m-3). Stations belonging to the low aCDOM(443) cluster (Cluster 1) are presented as circles, whereas stars represent stations grouped in the high aCDOM(443) cluster (Cluster 2; Fig 5). (b) Chl-a estimated by empirical algorithms (mg m-3; indicated by different symbols) versus in situ Chl-a (mg m-3), with aCDOM(443) (m-1) as colorbar. (c) In situ Chl-a (mg m-3) versus maximum band ratio [MBR; Rrs(443>490>510/555)]. (d) Chl-a estimated by empirical algorithms relative error (%) versus the ratio between aCDOM(443) and atw(443).

https://doi.org/10.1371/journal.pone.0190838.g007

thumbnail
Table 4. Evaluation of empirical ocean color algorithms.

Regression statistics (including the bias and the mean absolute error–MAE) for retrieved Chl-a from in situ Rrs compared to direct measurements of Chl-a using the correspondent algorithms versus in situ measured parameters. r2 and slope were calculated using log-transformed data for each of the correspondent parameters.

https://doi.org/10.1371/journal.pone.0190838.t004

A study in the western Arctic obtained good Chl-a retrievals applied to in situ Rrs measurements from CDOM-dominated waters (where Chl-a does not covary with CDOM), when turbid waters [Rrs(676)>0.00042] were excluded [35]. This could be one of the reasons attributed to the poor performance of those algorithms in our study, given that all the sampling stations were classified as turbid. This is supported by the fact that the most overestimated Chl-a retrievals were especially related to the high CDOM cluster (see Fig 7). When looking only at the stations grouped in Cluster I, i.e. with lower aCDOM(443) (although still with relatively high turbidity), the retrievals were significantly improved (strong reduction in RMSE), primarily for OC4V6-mod and OC3M-mod (RMSE = 0.08 and 0.10, respectively), and positively correlated to in situ Chl-a (Table 5). Estimates were less variable (r2>0.62) compared to previously reported for the western Arctic [35]. Since there are no specific empirical algorithms to derive CDOM in the Arctic, we have applied to our data set the algorithm developed by Belanger et al. [39] for the Western Arctic to obtain the ratio of aCDOM to adg in optically complex waters. However, we also did not obtain robust results using this method.

thumbnail
Table 5. Comparison of Chl-a retrieved from empirical ocean color algorithms versus direct measurements of Chl-a, for low aCDOM(443) sites.

Same as Table IV but for the low aCDOM(443) stations.

https://doi.org/10.1371/journal.pone.0190838.t005

Besides, the SAA GIOP was applied to retrieve Chl-a, as well as aph(443) and adg(443). This provided robust estimates (Fig 8) for the entire sampling area (Fig 8 and Table 6). Such an improvement probably is caused by that GIOP, like other SAAs, does not assume Chl-a and CDOM absorption as covariant. Even better estimates from GIOP were obtained for aph(443) (Fig 8 and Table 6). This is probably due to the fact that GIOP uses the spectral shape of Chl-a-specific absorption coefficient from Bricaud et al. [68] as basis vector. As pointed out before (Fig 3), our dataset exhibited similar spectral shape for aph(λ) and correlations between Chl-a and aph(443) as observed in that study. Moreover, the performance of GIOP to retrieve aph(443) in our study (Table 6) was much better than recently observed in the western Arctic (r2 = 0.85; Slope = 1.18; RMSE = 0.20) [33]. With regards to Chl-a, that same study reported fairly similar results (r2 = 0.72; Slope = 0.73; overestimation of Chl-a), however with lower errors (RMSE = 0.24) in comparison to our results (RMSE = 0.40, see Table 6).

thumbnail
Fig 8. Evaluation of semi-analytical algorithms.

Modeled geophysical parameters calculated from in situ Rrs versus in situ measured parameters: aph(443) (a); Chl-a (b); adg(443) (c) and aCDOM(443) (d). Red points refer to the GIOP [60,61] retrievals, whereas blue points to the retrievals from the GSM model adapted to the Arctic [19].

https://doi.org/10.1371/journal.pone.0190838.g008

thumbnail
Table 6. Evaluation of the semi-analytical ocean color algorithms.

Regression statistics for modeled geophysical parameters calculated used in situ Rrs versus in situ measured parameters. r2 and Slope were calculated using log-transformed data for each of the correspondent parameters.

https://doi.org/10.1371/journal.pone.0190838.t006

Finally, adg(443) was also retrieved in this study using GIOP and GSM-Matsuoka. Here we assume that adg(443) is a direct estimate from of aCDOM(443), given the least contribution of aNAP(443) to total non-water absorption (generally <1%) in comparison to aCDOM(443) (Fig 3). In general, GIOP and GSM-Matsuoka retrievals were very similar and the retrieved adg(443) compared very well to direct measurements of adg(443). Most of the data points located are within the 50% error intervals (Fig 8 and Table 6). Although with similar error (RMSE = 0.08), GIOP seems to provide more robust correlation (slope = 1.05) and less variable (r2 = 0.91) adg(443) estimates for the western Arctic, as observed for aph(443) [33]. As being highly correlated to adg(443), aCDOM(443) retrieved with GSM-Matsuoka provided very similar statistics related to adg(443) using that same model, thus resembling the observed with our in situ sampling (see Table 1). Compared to our study, better performance for retrieving aCDOM(443) using GSM-Matsuoka (r2 = 0.87; Slope = 0.97; RMSE = 0.07), was reported for the western Arctic using a much larger dataset [21]. Although the aCDOM(443) GSM-Matsuoka retrievals were not as good in our study, results show that SAAs in general have a high potential for obtaining reliable Chl-a estimates than empirical algorithms in high CDOM containing waters, besides the advantages of also providing other reliable retrievals such as adg(λ), aph(λ) and aCDOM(λ). Thus, products from SAAs are more suitable for application to biogeochemical studies in the Arctic Ocean, although improvement of the current algorithms is still requested, given the persistence of embedded errors to those retrievals, as demonstrated in this study.

4. Summary and outlook

Based on a quasi-synoptic sampling strategy over the surface Central-Eastern Arctic Ocean we reiterate the dominance of CDOM related to total non-water absorption through the entire region under study. As CDOM and DOC are strongly correlated in the Arctic Ocean [6,20,53], one can assume CDOM as a very reliable proxy for retrieving carbon concentrations in that basin. This can provide additional insight into the Arctic biogeochemical cycles. Our results show that aCDOM(443) and aph(443), together with temperature and salinity, are useful parameters for distinguishing hydrographic regimes within the Central Arctic. Despite the reduced number of sampling sites, hyperspectral AOPs retrieved from under water radiometric measurements were able to reproduce the major bio-optical features at the surface by differentiating between sites with low and high CDOM. As demonstrated for the Eastern Atlantic [54], bio-optical provinces efficiently reflected the ecosystem variability/biogeography proposed by Longhurst [92] and thus are a valuable tool for biogeochemical modeling. However, currently, practically the entire Arctic Ocean is still classified as a unique ecosystem unit, despite the reports of clear geographic patterns in different aspects [4,24,29,93,94]. Moreover, to our understanding, no study has proposed such a sub-division of the Boreal Polar Province (BPLR) into bio-optical units. Based on our findings we, therefore, propose a geographical characterization of the sampling regions into bio-optical provinces, which reflect hydrographic characteristics of the region with regard to the non-water absorption: Laptev Sea Shelf, Laptev Sea, Central Arctic/Transpolar Drift, Beaufort Gyre, and Eurasian/Nansen Basin. Moreover, it becomes clear that the characterization of provinces, in particular, in the highly seasonal variable Arctic Ocean, cannot hold true for every season and every year. Thus, although here we present a first, simple bio-optical classification, we recall that such variability has been observed along the Arctic and integrative biogeochemical studies would benefit from the advances in Arctic Ecosystem monitoring and management by improving the delimitation of such geographic units. Future perspectives using automated platforms (e.g., floats, ITP, gliders) with bio-optical (e.g., Chl-a and DOM fluorescence, and hyperspectral radiometry) and salinity sensors will allow to monitor the spatial and temporal variability within those biogeographic provinces.

The evaluation of empirical ocean color algorithms (including the regionally tuned ones) applied to our in situ Rrs(λ) measurements showed that those algorithms are inappropriate to estimate Chl-a in the Central-Eastern Arctic Ocean, exhibiting an overall inverse correlation with in situ Chl-a and a positive correlation with aCDOM(443). This reinforces the existence of a persuasive positive bias by CDOM absorption on empirical Chl-a estimates for the Arctic Ocean. The semi-analytical ocean color algorithm GIOP, on the other hand, retrieved reliable and less variable Chl-a estimates related to the empirical algorithms, as well as very good estimates for aph(443) and considerably well estimates of adg(443), as also reported to the western Arctic [33]. Fairly similar retrievals were obtained within the GSM model with the modifications for the Arctic Ocean [19] for adg(443). The better performance by SAAs is mainly attributed to the fact that these algorithms do not consider Chl-a and CDOM as covariant.

Finally, with the ongoing pressure of climate change over the Arctic environment, a better understanding on the dynamics of carbon stocks has been sought. Ocean color remote sensing appears to be a key tool on improving both the spatial and temporal monitoring of these stocks. However, accurate ocean color retrievals are required to get to a real estimate of stocks and processes involving organic and inorganic carbon in the Arctic Ocean. Thus, we recall that additional spectral bands would improve the performance of ocean color algorithms, as demonstrated for the GIOP in the western Arctic [33]. In addition, the coverage of ocean color remote sensing data in the Polar Regions needs to be increased (for our cruise we obtained no match-ups with satellite data) by investing in developing efficient atmospheric correction for adjacency effects and low illumination conditions.

Acknowledgments

We thank the captain and crew of the R/V Polarstern for their support during the ARKXXVI/3 cruise and the chief scientist Dr. Ursula Schauer. We acknowledge Prof. Dr. Antje Boetius, Dr. Alexandra Cherkasheva, Dr. Estelle Kilias and Olivia Serdeczny for RAMSES radiometric operation, water sampling and filtration during the cruise and Dr. Eva-Maria Nöthig, Christiane Lorenzen and Sonja Wiegmann (AWI) for their assistance in the preparation of the cruise. We further thank Sonja Wiegmann and Sandra Murawski for measurements of particulate absorption and pigment samples, respectively. Funding for the cruise was supplied partly by the Helmholtz Innovation Fund via the project Phytooptics (VH-NG-300). R.G.-A. was supported by a PhD fellowship from the Coordination for the Improvement of Higher Level Personnel (CAPES-Brazil, Grant 12362/12-3) in collaboration with the German Academic Exchange Service (DAAD).

References

  1. 1. Shiklomanov I. World fresh water resources. In: Gleick PH, editor. Water in Crisis: A guide to the World’s fresh water resources. New York: Oxford University Press; 1993. pp. 13–24.
  2. 2. Alling V, Sanchez-Garcia L, Porcelli D, Pugach S, Vonk JE, van Dongen B, et al. Nonconservative behavior of dissolved organic carbon across the Laptev and East Siberian seas. Global Biogeochem Cycles. 2010;24.
  3. 3. Fichot CG, Kaiser K, Hooker SB, Amon RMW, Babin M, Bélanger S, et al. Pan-Arctic distributions of continental runoff in the Arctic Ocean. Sci Rep. 2013;3: 1053. pmid:23316278
  4. 4. Lalande C, Nöthig E-M, Somavilla R, Bauerfeind E, Shevchenko V, Okolodkov Y. Variability in under-ice export fluxes of biogenic matter in the Arctic Ocean. Global Biogeochem Cycles. 2014;28: 571–583.
  5. 5. Sánchez-García L, Alling V, Pugach S, Vonk J, van Dongen B, Humborg C, et al. Inventories and behavior of particulate organic carbon in the Laptev and East Siberian seas. Global Biogeochem Cycles. 2011;25: n/a-n/a.
  6. 6. Walker SA, Amon RMW, Stedmon CA. Variations in high-latitude riverine fluorescent dissolved organic matter: A comparison of large Arctic rivers. J Geophys Res Biogeosciences. 2013;118: 1689–1702.
  7. 7. Schuur EAG, Abbott BW, Bowden WB, Brovkin V, Camill P, Canadell JG, et al. Expert assessment of vulnerability of permafrost carbon to climate change. Clim Change. 2013;119: 359–374.
  8. 8. Rabe B, Karcher M, Kauker F, Schauer U, Toole JM, Krishfield RA, et al. Arctic Ocean basin liquid freshwater storage trend 1992–2012. Geophys Res Lett. 2014;41: 961–968.
  9. 9. Rabe B, Dodd PA, Hansen E, Falck E, Schauer U, Mackensen A, et al. Liquid export of Arctic freshwater components through the Fram Strait 1998&ndash;2011. Ocean Sci. 2013;9: 91–109.
  10. 10. Serreze MC, Holland MM, Stroeve J. Perspectives on the Arctic’s Shrinking Sea-Ice Cover. Science (80-). 2007;315: 1533–1536. pmid:17363664
  11. 11. Stroeve JC, Kattsov V, Barrett A, Serreze M, Pavlova T, Holland M, et al. Trends in Arctic sea ice extent from CMIP5, CMIP3 and observations. Geophys Res Lett. 2012;39: 1–7.
  12. 12. Polyakov I V., Walsh JE, Kwok R. Recent changes of Arctic multiyear sea ice coverage and the likely causes. Bull Am Meteorol Soc. 2012;93: 145–151.
  13. 13. Nicolaus M, Katlein C, Maslanik J, Hendricks S. Changes in Arctic sea ice result in increasing light transmittance and absorption. Geophys Res Lett. 2012;39: n/a-n/a.
  14. 14. Li WKW, McLaughlin FA, Lovejoy C, Carmack EC. Smallest Algae Thrive As the Arctic Ocean Freshens. Science (80-). 2009;326: 539–539. pmid:19900890
  15. 15. Arrigo KR, van Dijken GL. Secular trends in Arctic Ocean net primary production. J Geophys Res. 2011;116: C09011.
  16. 16. Pabi S, van Dijken GL, Arrigo KR. Primary production in the Arctic Ocean, 1998–2006. J Geophys Res. 2008;113: C08005.
  17. 17. Hansell DA, Kadko D, Bates NR. Degradation of Terrigenous Dissolved Organic Carbon in the Western Arctic Ocean. Science (80-). 2004;304: 858–861. pmid:15131302
  18. 18. Bélanger S, Xie H, Krotkov N, Larouche P, Vincent WF, Babin M. Photomineralization of terrigenous dissolved organic matter in Arctic coastal waters from 1979 to 2003: Interannual variability and implications of climate change. Global Biogeochem Cycles. 2006;20: n/a-n/a.
  19. 19. Matsuoka A, Hooker SB, Bricaud A, Gentili B, Babin M. Estimating absorption coefficients of colored dissolved organic matter (CDOM) using a semi-analytical algorithm for southern Beaufort Sea waters: application to deriving concentrations of dissolved organic carbon from space. Biogeosciences. 2013;10: 917–927.
  20. 20. Stedmon CA, Amon RMW, Rinehart AJ, Walker SA. The supply and characteristics of colored dissolved organic matter (CDOM) in the Arctic Ocean: Pan Arctic trends and differences. Mar Chem. Elsevier B.V.; 2011;124: 108–118.
  21. 21. Matsuoka A, Babin M, Doxaran D, Hooker SB, Mitchell BG, Bélanger S, et al. A synthesis of light absorption properties of the Arctic Ocean: application to semianalytical estimates of dissolved organic carbon concentrations from space. Biogeosciences. 2014;11: 3131–3147.
  22. 22. Gonçalves-Araujo R, Granskog MA, Bracher A, Azetsu-Scott K, Dodd PA, Stedmon CA. Using fluorescent dissolved organic matter to trace and distinguish the origin of Arctic surface waters. Sci Rep. 2016;6: 33978. pmid:27667721
  23. 23. Granskog MA, Pavlov AK, Sagan S, Kowalczuk P, Raczkowska A, Stedmon CA. Effect of sea-ice melt on inherent optical properties and vertical distribution of solar radiant heating in Arctic surface waters. J Geophys Res Ocean. 2015;120: 7028–7039.
  24. 24. Fernández-Méndez M, Katlein C, Rabe B, Nicolaus M, Peeken I, Bakker K, et al. Photosynthetic production in the central Arctic Ocean during the record sea-ice minimum in 2012. Biogeosciences. 2015;12: 3525–3549.
  25. 25. Codispoti LA, Kelly V, Thessen A, Matrai P, Suttles S, Hill V, et al. Synthesis of primary production in the Arctic Ocean: III. Nitrate and phosphate based estimates of net community production. Prog Oceanogr. 2013;110: 126–150.
  26. 26. Ulfsbo A, Cassar N, Korhonen M, van Heuven S, Hoppema M, Kattner G, et al. Late summer net community production in the central Arctic Ocean using multiple approaches. Global Biogeochem Cycles. 2014;28: 1129–1148.
  27. 27. Okolodkov YB, Dodge JD. Biodiversity and biogeography of planktonic dinoflagellates in the Arctic Ocean. J Exp Mar Bio Ecol. 1996;202: 19–27.
  28. 28. Kilias E, Kattner G, Wolf C, Frickenhaus S, Metfies K. A molecular survey of protist diversity through the central Arctic Ocean. Polar Biol. 2014;37: 1271–1287.
  29. 29. Metfies K, Von Appen WJ, Kilias E, Nicolaus A, Nöthig EM. Biogeography and photosynthetic biomass of arctic marine pico-eukaroytes during summer of the record sea ice minimum 2012. PLoS One. 2016;11: 1–20. pmid:26895333
  30. 30. Stedmon CA, Granskog MA, Dodd PA. An approach to estimate the freshwater contribution from glacial melt and precipitation in East Greenland shelf waters using colored dissolved organic matter (CDOM). J Geophys Res Ocean. 2015;120: 1107–1117.
  31. 31. Matsuoka A, Babin M, Devred EC. A new algorithm for discriminating water sources from space: A case study for the southern Beaufort Sea using MODIS ocean color and SMOS salinity data. Remote Sens Environ. Elsevier Inc.; 2016;184: 124–138.
  32. 32. Ramírez-Pérez M, Gonçalves-Araujo R, Wiegmann S, Torrecilla E, Bardaji R, Röttgers R, et al. Towards cost-effective operational monitoring systems for complex waters: analyzing small-scale coastal processes with optical transmissometry. PLoS One. 2017;12: e0170706. pmid:28107539
  33. 33. Chaves JE, Werdell PJ, Proctor CW, Neeley AR, Freeman SA, Thomas CS, et al. Assessment of ocean color data records from MODIS-Aqua in the western Arctic Ocean. Deep Sea Res Part II Top Stud Oceanogr. Elsevier; 2015;118: 32–43.
  34. 34. Ben Mustapha S, Bélanger S, Larouche P. Evaluation of ocean color algorithms in the southeastern Beaufort Sea, Canadian Arctic: New parameterization using SeaWiFS, MODIS, and MERIS spectral bands. Can J Remote Sens. 2012;38: 535–556.
  35. 35. Matsuoka A, Huot Y, Shimada K, Saitoh S-I, Babin M. Bio-optical characteristics of the western Arctic Ocean: implications for ocean color algorithms. Can J Remote Sens. 2007;33: 503–518.
  36. 36. Reynolds R, Matsuoka A, Hirawake T, Mitchell BG. Ocean Colour algorithms and bio-optical relationships for polar seas, IOCCG Report 16. 2015.
  37. 37. IOCCG. IOCCG Report number 16: Ocean Colour Remote Sensing in Polar Seas. 2015.
  38. 38. Cota G, Wang J, Comiso JC. Transformation of global satellite chlorophyll retrievals with a regionally tuned algorithm. Remote Sens Environ. 2004;90: 373–377.
  39. 39. Bélanger S, Babin M, Larouche P. An empirical ocean color algorithm for estimating the contribution of chromophoric dissolved organic matter to total light absorption in optically complex waters. J Geophys Res. 2008;113: C04027.
  40. 40. Krishfield R, Toole J, Proshutinsky A, Timmermans ML. Automated ice-tethered profilers for seawater observations under pack ice in all seasons. J Atmos Ocean Technol. 2008;25: 2091–2105.
  41. 41. Laney SR, Krishfield RA, Toole JM, Hammar TR, Ashjian CJ, Timmermans ML. Assessing algal biomass and bio-optical distributions in perennially ice-covered polar ocean ecosystems. Polar Sci. 2014;8: 73–85.
  42. 42. Toole J, Krishfield R, Timmermans M-L, Proshutinsky A. The Ice-Tethered Profiler: Argo of the Arctic. Oceanography. 2011;24: 126–135.
  43. 43. Schauer U, Rabe B, Wisotzki A. Physical oceanography during POLARSTERN cruise ARK-XXVI/3 (TransArc). PANGAEA. 2012;
  44. 44. Schlitzer R. Ocean Data View [Internet]. 2015. Available: http://odv.awi.de
  45. 45. Rudels B. Arctic Ocean Circulation. In: Steele JH, Thorpe SA, Turekian KK, editors. Ocean Currents: A Derivative of the Encyclopedia of Ocean Sciences. 2009. pp. 211–225. Available: https://books.google.com/books?id=FYSCUH235E8C&pgis=1
  46. 46. Simis SGH, Tijdens M, Hoogveld HL, Gons HJ. Optical changes associated with cyanobacterial bloom termination by viral lysis. J Plankton Res. 2005;27: 937–949.
  47. 47. Röttgers R, Gehnke S. Measurement of light absorption by aquatic particles: improvement of the quantitative filter technique by use of an integrating sphere approach. Appl Opt. 2012;51: 1336–1351. pmid:22441480
  48. 48. Tassan S, Ferrari GM. An alternative approach to absorption measurements of aquatic particles retained on filters. Limnol Oceanogr. 1995;40: 1358–1368.
  49. 49. Ferrari GM, Tassan S. A method using chemical oxidation to remove light absorption by phytoplankton pigments. J Phycol. 1999;35: 1090–1098.
  50. 50. Heim B, Abramova E, Doerffer R, Günther F, Hölemann J, Kraberg A, et al. Ocean colour remote sensing in the southern Laptev Sea: evaluation and applications. Biogeosciences. 2014;11: 4191–4210.
  51. 51. Siegel DA, Maritorena S, Nelson NB, Behrenfeld MJ, McClain CR. Colored dissolved organic matter and its influence on the satellite-based characterization of the ocean biosphere. Geophys Res Lett. 2005;32: L20605.
  52. 52. Prieur L, Sathyendranath S. An optical classification of coastal and oceanic waters based on the specific spectral absorption curves of phytoplankton pigments, dissolved organic matter, and other particulate materials1. Limnol Oceanogr. 1981;26: 671–689.
  53. 53. Gonçalves-Araujo R, Stedmon CA, Heim B, Dubinenkov I, Kraberg A, Moiseev D, et al. From Fresh to Marine Waters: Characterization and Fate of Dissolved Organic Matter in the Lena River Delta Region, Siberia. Front Mar Sci. 2015;2: 108.
  54. 54. Taylor BB, Torrecilla E, Bernhardt A, Taylor MH, Peeken I, Röttgers R, et al. Bio-optical provinces in the eastern Atlantic Ocean and their biogeographical relevance. Biogeosciences. 2011;8: 3609–3629.
  55. 55. Jeffrey SW, Mantoura RFC, Wright SW, editors. Phytoplankton pigments in Oceanography: Guideline to modern methods. UNESCO Publishing; 1997.
  56. 56. Mueller JL, Fargion GS, McClain CR. Ocean Optics Protocols for Satellite Ocean Color Sensor Validation, Revision 4, Volume III: Radiometric Measurements and Data Analysis Protocols. Greenbelt, Maryland: NASA/TM-2003-211621/Rev4-Vol. III, NASA Goddard Space Flight Center; 2003.
  57. 57. Stramski D, Reynolds RA, Babin M, Kaczmarek S, Lewis MR, Röttgers R, et al. Relationships between the surface concentration of particulate organic carbon and optical properties in the eastern South Pacific and eastern Atlantic Oceans. Biogeosciences. 2008;5: 171–201.
  58. 58. O’Reilly JE, Maritorena S, O’Brien MC, Siegel DA, Toole D, Menzies D, et al. SeaWiFS Postlaunch Calibration and Validation Analyses, Part 3. SeaWiFS Postlaunch Tech Rep Ser. 2000;11: 51pp.
  59. 59. O’Reilly JE, Maritorena S, Mitchell BG, Siegel DA, Carder KL, Garver SA, et al. Ocean color chlorophyll algorithms for SeaWiFS. J Geophys Res Ocean. 1998;103: 24937–24953.
  60. 60. Werdell PJ, Franz BA, Bailey SW, Feldman GC, Boss E, Brando VE, et al. Generalized ocean color inversion model for retrieving marine inherent optical properties. Appl Opt. 2013;52: 2019–2037. pmid:23545956
  61. 61. Werdell PJ, Franz BA, Lefler JT, Robinson WD, Boss E. Retrieving marine inherent optical properties from satellites using temperature and salinity-dependent backscattering by seawater. Opt Express. 2013;21: 32622. pmid:24514855
  62. 62. Garver SA, Siegel DA. Inherent optical property inversion of ocean color spectra and its biogeochemical interpretation: 1. Time series from the Sargasso Sea. J Geophys Res Ocean. 1997;102: 18607–18625.
  63. 63. Maritorena S, Siegel DA, Peterson AR. Optimization of a semianalytical ocean color model for global-scale applications. Appl Opt. 2002;41: 2705. pmid:12027157
  64. 64. Clarke KR, Warwick RM. Changes in Marine Communities: An approach to Statistical Analysis and Interpretation. Plymouth, UK: Natural Environmental Research Council; 1994.
  65. 65. Torrecilla E, Stramski D, Reynolds RA, Millán-Núñez E, Piera J. Cluster analysis of hyperspectral optical data for discriminating phytoplankton pigment assemblages in the open ocean. Remote Sens Environ. Elsevier Inc.; 2011;115: 2578–2593.
  66. 66. Campbell JW. The lognormal distribution as a model for bio-optical variability in the sea. J Geophys Res. 1995;100: 13237.
  67. 67. Bricaud A, Babin M, Morel A, Claustre H. Variability in the chlorophyll-specific absoprtion coefficients of natural phytoplankton: Analysis and parametrization. J Geophys Res. 1995;100: 13321–13332.
  68. 68. Bricaud A, Morel A, Babin M, Allali K, Claustre H. Variations of light absorption by suspended particles with chlorophyll a concentration in oceanic (case 1) waters: Analysis and implications for bio-optical models. J Geophys Res Ocean. 1998;103: 31033–31044.
  69. 69. Bauch D, Cherniavskaia E, Timokhov L. Shelf basin exchange along the Siberian continental margin: Modification of Atlantic Water and Lower Halocline Water. Deep Sea Res Part I Oceanogr Res Pap. Elsevier; 2016;115: 188–198.
  70. 70. Jones EP, Nelson DM, Treguer P. Chemical Oceanography. In: Smith WO Jr, editor. Polar Oceanography Part B: Chemistry, Biology, and Geology. Elsevier Inc.; 1990. pp. 407–476.
  71. 71. Harms IH, Karcher MJ. Modeling the seasonal variability of hydrography and circulation in the Kara Sea. J Geophys Res. 1999;104: 13431.
  72. 72. Janout MA, Aksenov Y, Hölemann JA, Rabe B, Schauer U, Polyakov I V., et al. Kara Sea freshwater transport through Vilkitsky Strait: Variability, forcing, and further pathways toward the western Arctic Ocean from a model and observations. J Geophys Res Ocean. 2015;120: 4925–4944.
  73. 73. Bauch D, Dmitrenko I, Kirillov S, Wegner C, Hölemann J, Pivovarov S, et al. Eurasian Arctic shelf hydrography: Exchange and residence time of southern Laptev Sea waters. Cont Shelf Res. 2009;29: 1815–1820.
  74. 74. Stramska M, Stramski D, Hapter R, Kaczmarek S, Stoń J. Bio-optical relationships and ocean color algorithms for the north polar region of the Atlantic. J Geophys Res. 2003;108: 3143.
  75. 75. Stramska M, Stramski D, Kaczmarek S, Allison DB, Schwarz J. Seasonal and regional differentiation of bio-optical properties within the north polar Atlantic. J Geophys Res. 2006;111: C08003.
  76. 76. Kowalczuk P, Meler J, Kauko HM, Pavlov AK, Zabłocka M, Peeken I, et al. Bio-optical properties of Arctic drift ice and surface waters north of Svalbard from winter to spring. J Geophys Res Ocean. 2017;122: 4634–4660.
  77. 77. Lund-Hansen LC, Markager S, Hancke K, Stratmann T, Rysgaard S, Ramløv H, et al. Effects of sea-ice light attenuation and CDOM absorption in the water below the Eurasian sector of central Arctic Ocean (>88°N). Polar Res. 2015;34: 1–12.
  78. 78. Granskog MA, Stedmon CA, Dodd PA, Amon RMW, Pavlov AK, de Steur L, et al. Characteristics of colored dissolved organic matter (CDOM) in the Arctic outflow in the Fram Strait: Assessing the changes and fate of terrigenous CDOM in the Arctic Ocean. J Geophys Res. 2012;117: C12021.
  79. 79. Pavlov AK, Granskog MA, Stedmon CA, Ivanov B V, Hudson SR, Falk-Petersen S. Contrasting optical properties of surface waters across the Fram Strait and its potential biological implications. J Mar Syst. Elsevier B.V.; 2015;143: 62–72.
  80. 80. Børsheim KY, Myklestad SM. Dynamics of DOC in the Norwegian Sea inferred from monthly profiles collected during 3 years at 66°N, 2°E. Deep Sea Res Part I Oceanogr Res Pap. 1997;44: 593–601.
  81. 81. Siegel DA, Maritorena S, Nelson NB, Hansell DA, Lorenzi-Kayser M. Global distribution and dynamics of colored dissolved and detrital organic materials. J Geophys Res. 2002;107: 3228.
  82. 82. Nelson NB, Siegel DA. The Global Distribution and Dynamics of Chromophoric Dissolved Organic Matter. Ann Rev Mar Sci. 2013;5: 447–476. pmid:22809178
  83. 83. Morel A, Prieur L. Analysis of variations in ocean color. Limnol Oceanogr. 1977;22: 709–722.
  84. 84. Antoine D, Babin M, Berthon J-F, Bricaud A, Gentili B, Loisel H, et al. Shedding Light on the Sea: André Morel’s Legacy to Optical Oceanography. Ann Rev Mar Sci. 2014;6: 1–21. pmid:24015899
  85. 85. Cota GF, Harrison WG, Platt T, Sathyendranath S, Stuart V. Bio-optical properties of the Labrador Sea. J Geophys Res. 2003;108: 3228.
  86. 86. Wang J, Cota GF, Ruble DA. Absorption and backscattering in the Beaufort and Chukchi Seas. J Geophys Res. 2005;110: C04014.
  87. 87. Babin M, Stramski D, Ferrari GM, Claustre H, Bricaud A, Obolensky G, et al. Variations in the light absorption coefficients of phytoplankton, nonalgal particles, and dissolved organic matter in coastal waters around Europe. J Geophys Res. 2003;108: 3211.
  88. 88. Semiletov IP, Shakhova NE, Pipko II, Pugach SP, Charkin AN, Dudarev O V., et al. Space-time dynamics of carbon and environmental parameters related to carbon dioxide emissions in the Buor-Khaya Bay and adjacent part of the Laptev Sea. Biogeosciences. 2013;10: 5977–5996.
  89. 89. Proshutinsky A, Krishfield R, Timmermans M-L, Toole J, Carmack E, McLaughlin F, et al. Beaufort Gyre freshwater reservoir: State and variability from observations. J Geophys Res. 2009;114: C00A10.
  90. 90. Bluhm BA, Kosobokova KN, Carmack EC. A tale of two basins: An integrated physical and biological perspective of the deep Arctic Ocean. Prog Oceanogr. Elsevier Ltd; 2015;139: 89–121.
  91. 91. Neukermans G, Reynolds RA, Stramski D. Optical classification and characterization of marine particle assemblages within the western Arctic Ocean. Limnol Oceanogr. 2016;61: 1472–1494.
  92. 92. Longhurst A. Ecological geography of the sea. 2nd ed. Elsevier Inc.; 2007.
  93. 93. Arrigo KR, Hill V, Bélanger S, Mitchell BG, Hirawake T, Babin M. Estimates of Net Primary Production from Space-based Measurements, IOCCG Report 16: Ocean Colour Remote Sensing in Polar Seas. 2015.
  94. 94. Rutgers van der Loeff MM, Cassar N, Nicolaus M, Rabe B, Stimac I. The influence of sea ice cover on air-sea gas exchange estimated with radon-222 profiles. J Geophys Res Ocean. 2014;119: 2735–2751.